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Abstract. We have systematically estimated the possible temperatures obtained from an analysis of recent 
data on pt distributions observed at RHIC experiments. Using the fact that observed pt distributions cannot 
be described by the original Hagedorn formula in the whole range of transverse momenta (in particular 
above 6 GeV/c), we propose a modified Hagedorn formula including temperature fluctuation. We show 
that by using it we can fit pt distributions in the whole range and can estimate consistently the relevant 
temperatures, including their fluctuations. 

PACS. 25.75.-q Relativistic heavy ion collisions - 12.40.Ee Statistical (extensive and non-extensive) models 
- 02.50.Ey Stochastic models 



1 Introduction 

One of the characteristic features in every high energy 
collision experiment is the production of large numbers 
of secondaries (mostly pions). From the very beginning of 
the history of the multiparticle production processes, it 
was realized that a possible way to treat them was to em- 
ploy some sort of statistical approach P . This idea found 
its most mature formulation in the statistical bootstrap 
model proposed by Hagedorn j2, , in which the exponential 
growth of the number of hadronic resonances with mass 
is one of the most fundamental issues 1^. The proposed 
formula is 



data on resonance production ^ . The other parameter is 
Po = l/(fc_BTo), with To explicitly governing the observed 
energy distribution and therefore identified with the tem- 
perature of the hadronizing system. In the foUowings we 
put fcs = 1. One of the aims in the study of multipar- 
ticle production processes is therefore the best possible 
estimation of this quantity. To this end we would like to 
investigate the measured transverse momentum (pt) dis- 
tributions integrated over longitudinal degrees of freedom. 
From Eq.(^, we have 



^0 = 



2'Kptdpt 



= C dmp{m)mtKi{mtPo) 



(3) 



fa 
dp^ 

In Eq.C 

by 



= C 



dmp{m) exp (^-^pf + P^ + m'^l^o^ ■ (1) 



p(m) denotes the density of resonances given 



pirn) 



exp(m/?g) 



where fin = ^/{kBTH), the inverse of the so called Hage- 
dorn temperature Th, is a parameter to be deduced from 
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where mt = \J p\ ^ m? is transverse mass and Ki is the 
Bessel function. 

However, as was recently demonstrated by us l^j, this 
simple formula can explain the RHIC data only in the 
limited range of transverse momenta, namely for < 6 
GeV/c. For larger values of pt data exhibit a power-like 
tail. There are many attempts to explain it using some 
kind of nonequilibrium approach like, for example, the 
flow or decay of resonances (see [HI for most recent re- 
view and further references); instead of trying to exclude 
them we would like to investigate the possibility that the 
observed nonexponential spectra could result from some 
form of equilibrium characteristic of nonextensive thermo- 
dynamics. In fact, as was shown in [2], using an approach 
based either on nonextensive statistics or on stochastic ap- 
proach one can successfully account for the whole range of 
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the observed transverse niomenta. The reason for this suc- 
cess is the fact that in both approaches the resuhant dis- 
tributions are intrinsically non- exponential, ranging from 
a power-law like form (cf. Eq.® below) to a gaussian 
in transverse rapidity [7] (which can be regarded as an- 
other implementation of the effective power-law distribu- 
tion)^. The fact that the proposed formulas can fit the 
whole range of pt is by itself very interesting and impor- 
tant observation as it shows that the power-law is present 
not only in very hard scale physics but that it reflects 
also a possible nontrivial property of hadronic matter in 
equilibrium (like, for example, Quark Gluon Plasma) j^i- 
Such properties are best seen in an approach using a 
nonextensive statistical model in which two parameters 
are used: the action of the heat bath is described now 
by the mean temperature Tq and by the nonextensivity 
parameter q, which can be identified with some specific 
intrinsic fluctuations of the temperature existing in the 
hadronizing system under consideration |10j . In the case 
when these fluctuations can be described by gamma dis- 
tribution one can write exact formulas TIP telling us that 

(cf. dm) 



[l-(l-g)/3offo]^ 

where 

1 / a \ " 



e-^"«Mf3)d(3, (4) 



r(a) V/3o 



9-1 



(5) 



In general, one refers to the concept of so called super- 
statistics introduced in . In our previous work , RHIC 
data were described by the following distribution with 



rnl: 



(Pa 



2'Kptdpt 



C 



dpi 



1-g 
To 



Pt 



ml 



(6) 



As is seen in , this formula leads to very good agreement 
with all RHIC data pill^lll4| . 

It is important to notice that Eq. has essentially 
the same form as the formula proposed long time ago and 
used with success in many QCD-inspired power- law fits to 
experimental data jl5lll6lll7j (recently used also by RHIC 
collaborations [TH|): 



1 + ^ 



exp 



(-^p,) forp, 
for Pt — !■ oo. 



(7) 



However, one has also to realize the important difference 
in physical pictures leading to Eq. (|HJ and Eq. 0. The 
underlying physical picture in Eq. Q is that the small 
Pt region is governed by soft physics described by some 
unknown unperturbative theory or model, and the large 

^ See also where flow effect is included and relation be- 
tween Gaussian-like distribution in transverse rapidity and 
power law behavior in pt is discussed. 



Pt region is governed by hard physics believed to be de- 
scribed by perturbative QCD. Contrary to it, the nonex- 
tensive formula Eq.®, which is valid in the whole range 
of Pt , does not claim to originate from any particular the- 
ory. It merely offers the kind of general unifying principle, 
namely the existence of some kind of complicated equi- 
librium involving all scales of pt, which is described by 
two parameters, Tq and q: the temperature Tq describing 
its mean properties and the parameter q describing ac- 
tion of the possible nontrivial long range effects believed 
to be caused by fluctuations but essentially also by some 
correlations or long memory effects ^• 



2 Calculations and results 

In this paper, we would like to compare results of an analy- 
sis of spectra measured at RHIC experiments |12II13II14| 
performed by using three approaches: the original Hage- 
dorn model, Eq. Q, the QCD-inspired power-like formula, 
Eq. |(7J), and the modified Hagedorn formula including 
temperature fluctuation given by : 



C J dy coshy J dmp{m)mt 



2-Kptdpt 

•[l-/3o(l-g)mtCoshy]T^ . (8) 
It can be written also in the form of series (a = \/{q—l)): 



St, 



1 



dmp{m) 



(1 + Pomt/a)" 



(9) 



■E 

fe=0 



r{k + 3/2)r(a + 1 + fc) /I - Pomt/a 



r{a + k + l/2)r{k +1) V 1 + Pomt/a 



or, accounting for the smallness of g — 1 encountered in our 
fits and of the fact that we are interested only in midra- 
pidity region (i.e., for small y) one can write it also as^ 



Stot — C dy coshy / dmp{m)mt 



1 



■ exp 



-^o"it coshy +-((?- l)/3o"^t 



(10) 



Equations © and ((117)) are used to check the numerical 
integration of Eq. © . 

At first the STAR data ^21 were analyzed using of 
Eq. JSJ (which corresponds to g = 1 in Eq. ©) and modi- 
fied Hagedorn formula, Eq. ((SJ. Results are shown in Fig.^ 

^ The origin of such fluctuations and/or correlations must be 
most probably traced back to the nonperturbative QCD, cf., 
for example, |2()| . 

^ In our case, because we are integrating over the whole mass 
spectrum pirn) in Hagedorn formula, we cannot simply expand 
in (g — 1) and keep only linear term as it is done on such 
occasions in the literature, cf., for example |21| . because for 
large masses m the series becomes divergent. 
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Fig. 1. Analysis of STAR data 12 by using usual Hagedorn formula (Eq.Q, left panel) and its nonextensive generalization 
(Eq. JSJ, right panel). 



Table 1. Parameters of our analysis presented in Fig. Q (left panel) by the use of Eq. Q which corresponds to g — 1 = in 
Eq.JHJ. Those for right panel with g — 1 7^ 1 in Eq.® can be found in Table 4. Other parameters are mo = 0.5 GeV (fixed), 
STh = 0.0001-0.002 and 5Tq = 0.0001-0.002. Notice that very large values of x^are obtained for fits with g - 1 = 0. 



C. C. (%) 


C 


Th (GeV) 


To (GeV) 


xVn.d.f. 


0-5 


816± 15 


0.086 


0.085 


532/32 


20-30 


382±7 


0.077 


0.076 


249/32 


60-80 


106± 2 


0.037 


0.037 


308/32 




11.35 11.4 11.45 11.5 6.64 6.66 6.68 6.7 6.72 6.74 6.76 6.78 

j8 [GeV-i] P [GeV"'] 



Fig. 2. Temperature fluctuation in STAR data (for C.C. = — 5 %) ^2] are analyzed by the use of Eq. Q with g — 1 = (left 
panel) and Eq. with g — 1 7^ (right panel). 



and Tabled The corresponding results for BRAHMS and 
PHENIX data JJ. 14 are very similar. As can be seen in 
Fig. ^ whereas distributions in the small pt region can 
be explained by simple formula (jSJ, data including the 
larger pt region can only be explained by using the modi- 
fied Hagedorn formula, Eq. JHl (or Eqs. © and ifTUI) ). The 
nonzero values of \q — 1| are then interpreted as an indi- 
cation of the sizeable temperature fluctuations existing in 
the hadronizing system 

In fact RHIC data allow us to investigate the fluctua- 
tion of temperature in more detail, cf. Fig.|21 The central- 
ity cut C.C. = - 5 % region of STAR data 52 was fitted 



by using, respectively, Eq. (with q — 1 — Q, left panel) 
and Eq. |(SJ) (with q as in Table 4, right panel). Fit was per- 
formed by fixing all parameters in Eqs. @ and (O except 
/So; (3{— /3o) is then calculated for each of 35 data points 
and it is assumed that the reciprocal of the each error 
bar calculated by the fitting program MINUIT is propor- 
tional to the corresponding probability of this value of /?, 
P(/3). In this way a probability distribution for (3 is ob- 
tained and presented in the form of a histogram in Fig|21 
The histogram in each panel is then fitted to the Gamma 
distribution with a = 55000 (shown as solid curves). The 
mean value (/?) is also shown in each panel in FigEl As can 
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Fig. 3. Example of tlie visualization of results presented in Tables 2 and 4 using STAR |12| results. Left panels sliow in detail 
contribution of different medianisms represented by Stot Elq. l|HJ , by 5o Eq. @ and by tfieir difference denoted by Sq — Stot — So- 
Rigfit panels sfiow results on original Hagedorn model, Eq.Q by daslied line, the modified the Hagedorn formula including 
temperature fluctuation, Eq.||HJ, by solid line, and Eq.JTJ by dotted fine. 



be seen, a good fit can be obtained only when Eq. © is 
used and in this case the resultant distribution of temper- 
atures is very narrow. This result suggests that accounting 
for intrinsic fluctuations considerably narrows the distri- 
bution of temperatures (actually its reverse, (3 = 1/T) and 
minimizes what can be regarded as a kind of systematic 
error in deduction of /So from experimental data. Therefore 
it strongly suggests that the modified Hagedorn formula, 
Eq. should be used whenever possible. 

The results of our fits to RHIC data VTJTJ^ per- 
formed by using Eq. ® (as given by noncxtensivc statis- 
tical approach), Eq. O (representing the QCD-inspired 
power-law formula) and Eq. Q (given by the modified 



Hagedorn formula proposed by us here) are presented in, 
respectively. Tables 2, 3 and 4. The results for STAR 
data are then also shown in Fig|31 In particular, the left 
hand panels of Fig.|21demonstrate contribution of different 
mechanism represented, respectively by Stot and 5*0. It is 
clear that data for the larger pt region can be explained 
only by Stot , which can be attributed to the intrinsic pri- 
mordial temperature fluctuations in the hadronizing sys- 
tem. However, at present it is difficult to treat this as 
a possible signal of a Quark-Gluon Plasma. Notice that 
the temperature parameter Tq ~ l//3o in Table 2 and 4 
was estimated by the use of Eq. (jSl from the whole re- 
gion of transverse momenta, whereas Tq = Po/n^ which 
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corresponds to temperature in Eq. 0, shown in Table 3 
governs only the small pt region. RHIC data show that 
we always have To > To, i.e., that inclusion of fluctuations 
and long-range correlations present in the hadronizing sys- 
tem lowers the estimated value of its mean temperature. 
From Table 4, we can see that both temperatures, Th 
and To, estimated by the use of Eq. JHJ decrease as the 
centrality cut, C.C., increases (i.e., it can be argued that 
they increase with the volume of interaction; similar ef- 
fect concerning Th has been also found in ^TT). It should 
be emphasized that when one uses the modified Hagedorn 
formula, Eq. (|2J), then Th ^ Tq TOtt, i-e., estimated val- 
ues of Th and Tq are almost equal to m^r, which we regard 
as very reasonably result*. 

3 Summary 

We have presented a systematic analysis of RHIC data [121 
I13II14| on transverse momenta distributions, which allow, 
in principle, the deduction of the parameter believed to 
represent the temperature Tq of the hadronizing system. 
We have shown that in order to fit the whole range of pt 
one has to use a nonextensive approach, which accounts 
for temperature fluctuations present in the hadronizing 
system. This has been compared with approach using the 
old QCD-inspired power-like formulas introduced long time 
ago. We have demonstrated that gradual accounting for 
the intrinsic dynamical fluctuations in the hadronizing 
system by switching from Eq. © (as given by nonexten- 
sive statistical approach) to the modified Hagedorn for- 
mula including temperature fluctuation, Eq. (jHJ , substan- 
tially lowers the values of parameter q—1. This is because 
part of the fluctuations ascribed in Eq. © to g are ac- 
counted for by the resonance spectrum p{m) present in 
the Hagedorn formula. It also changes the temperature 
we are looking for. Therefore one has to be very careful 
when interpreting the temperature parameter obtained in 
such fits, especially when attempting to address any ques- 
tions concerning Quark Gluon Plasma production issues^ . 

Actually, analysis performed assuming both thermal and 
chemical equilibrium and including also baryons performed by 
GSI group '23' gives T = 170 MeV. In our case we are consid- 
ering only pious and get T ~ m^. This difference is important 
for the description of phase diagram and we plan to address it 
elsewhere. One more remark is in order here. The To parame- 
ters obtained by us are in a range of Tc — 170 MeV, the QCD 
crossover temperature. On the other hand, traditional expo- 
nential fits for the low pr part of pion spectra used to give 
T = 340 MeV, pointing to a transverse flow with a Doppler 
blue-shift factor of two. However, we do not claim that there 
is no transverse flow in RHIC experiments, we only show that 
nonextensive approach can mimic this effect as well. 

^ One should be aware of the fact that there is still an on- 
going discussion on the meaning of the temperature in nonex- 
tensive systems. However, the small values of the parameter 
q deduced from data allow us to argue that, to first approxi- 
mation. To can be regarded as the hadronizing temperature in 
such system. One must only remember that in general what we 
study here is not so much the state of equilibrium but rather 



If data with larger pt are available, we can further inves- 
tigate whether the modifled Hagedorn formula including 
temperature fluctuation is really applicable or not. 
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Table 2. Analysis of RHIC data |12ll3ll3] by means of nonextensive approach as given by Eq. Q. For comparison the results 
for PP collisions are also shown. 



CnU C, C, 






9 


^ /(n -^) 


/n H f 
X /n.u.i 


ijrv,/\niVio u-iu/o 




n 997-i-n nn^ 

\J.LL t ZIZU.UUO 


1 o'^Qzi+0 nn9fi 


9^1 d 


10 9 /99 


lU-zu /o 


t lOICOD 


n 91 7-1-0 nn^ 


1 0/1 ^Pi^o nn9Q 


99 


1 9 Q /9^ 


zu-^u /O 




n 908-1-0 oofi 


1 0^07^0 00'^'^ 


1 Q 7 

iy. ( 


1 9 S /9^ 
iz.o / zo 


zLO (S0% 




n 1 s'i-i-o 007 


1 ofi07+o nrizizi 




1 fi /9^ 


an eo% 


Q« 9-1-4. 9 


1 fi^i-l-O 00^ 


1 07fizl+0 0094- 


1 1 


9 7fi/9^ 




1 c;Qn-i-Qc;Q 
iOoUztoOy 


1 Qf^-I-O 019 


1 0/1fi1 -1-0 ClClf^Cl 


91 7 


c; no /9Q 
o.uu/ zy 


0-10 /o 




1 QQ-t-O 019 


1 0/179-1-0 00f^7 


91 9 


o rf; /9Q 

o.OD/ zy 


15-S0% 


760+1 80 


lSq+0 012 


1 0503+0 0058 


19.9 


5 50 /29 


30-60% 


384±96 


0.170±0.011 


1.0613+0.0055 


16.3 


2.60/29 


60-80% 


120±39 


0.144±0.012 


1.0728+0.0067 


13.7 


10.5/29 


80-92% 


59.2±32.0 


0.114±0.017 


1.0879+0.0106 


11.4 


8.99/29 


STAR 0-5% 


3980±186 


0.164±0.002 


1.0651+0.0009 


15.4 


172/32 


5-10% 


2900±148 


0.169±0.002 


1.0622+0.0011 


16.1 


64.5/32 


10-20% 


2340±114 


0.164±0.002 


1.0662+0.0011 


15.1 


66.4/32 


20-30% 


1630±81 


0.162±0.002 


1.0684+0.0011 


14.6 


40.7/32 


30-40% 


1170±61 


0.158±0.002 


1.0709+0.0011 


14.1 


38.9/32 


40-60% 


739±39 


0.146±0.002 


1.0772+0.0010 


13.0 


14.7/32 


60-80% 


328±19 


0.130±0.002 


1.0850+0.0011 


11.8 


9.39/32 


PP (nsd) 


49.9±5.5 


0.111±0.003 


1.0894+0.0014 


11.2 


10.1/29 



Table 3. Analysis of RHIC data |12ll3lll] by means of the QC- inspired power- like formula Q. For comparison the results for 
PP collisions are also shown. 



CoU. & C.C. 


c 


n 


Po [GeV] 


r = po/n 


xVn.d.f 


BRAHMS 0-10% 


353+19 


32.2+3.3 


8.89+1.09 


0.276 


7.96/23 


10-20% 


260+15 


26.4+2.5 


7.05+0.83 


0.267 


14.7/23 


20-40% 


163+11 


22.8+2.1 


5.87+0.70 


0.257 


13.9/23 


40-60% 


83.7+7.5 


17.9+1.8 


4.13+0.56 


0.231 


11.7/23 


60-80% 


11.1+1.0 


12.8+0.5 


2.58+0.17 


0.202 


2.86/23 


PHENIX 0-5% 


536+398 


23.8+34.0 


5.54+31.26 


0.233 


4.69/29 


5-15% 


417+276 


23.0+23.7 


5.30+21.15 


0.231 


3.58/29 


15-30% 


260+149 


21.3+13.2 


4.84+11.01 


0.227 


5.54/29 


30-60% 


120+924 


16.6+91.8 


3.40+90.58 


0.205 


2.66/29 


60-80% 


32.1+38.2 


13.5+31.5 


2.38+30.46 


0.177 


10.3/29 


80-92% 


12.8+35.4 


10.8+31.6 


1.53+30.44 


0.142 


8.83/29 


STAR 0-5% 


1140+41 


15.4+0.3 


3.10+0.09 


0.201 


194/32 


5-10% 


843+33 


16.5+0.4 


3.44+0.12 


0.208 


68.6/32 


10-20% 


660+25 


15.3+0.3 


3.12+0.10 


0.203 


72.3/32 


20-30% 


457+18 


14.7+0.3 


2.94+0.09 


0.200 


42.8/32 


30-40% 


319+13 


14.1+0.3 


2.77+0.09 


0.196 


38.1/32 


40-60% 


190+8 


12.6+0.2 


2.30+0.07 


0.182 


13.9/32 


60-80% 


75.4+3.3 


11.3+0.2 


1.84+0.06 


0.163 


7.18/32 


PP (nsd) 


10.8+0.9 


10.4+0.2 


1.42+0.06 


0.136 


11.6/29 
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Table 4. Analysis of RHIC data [1211811^ by means of nonextensive modification of the Hagedorn formula as given by Eq. Q. 
Maximum m is fixed at 70 GeV (therefore in @ one always has (1 — Pomt/a) > 0). Numbers of divisions for y and m in 
computations are given in the last column. For comparison the results for pp collisions are also shown. 



CoU. & C.C. 


c 


q-l 




T„ [GeV] 


To [GeV] 


X7n.d.f. 


no. of div. 


BRAHMS 0-10% 


156±3 


0.00 




0.192±0.000 


0.178±0.000 


15.4/22 


6x6 


10-20% 


106±5 


(4.76±0.62)xl0" 


-4 


0.206±0.007 


0.187±0.005 


13.1/22 


6x7 


20-40% 


67.7±4.9 


(8.49±21.05)xl0 


-5 


0.177±0.013 


0.166±0.010 


11.6/22 


5x3 


40-60% 


32.5±2.9 


(2.57±0.63)xl0" 


-4 


0.168±0.010 


0.157±0.008 


9.54/22 


5x4 


60-80% 


5.00±0.14 


(8.12±0.45)xl0~ 


-5 


0.124±0.000 


0.120±0.000 


3.19/22 


6x6 


PHENIX 0-5% 


226±56 


(1.21±2.31)xl0~ 


^4 


0.16±0.02 


0.152±0.019 


4.98/29 


5x6 


5-15% 


157±34 


(4.01±0.02)xl0~ 


-4 


0.183±0.023 


0.167±0.017 


3.32/29 


6x5 


15-30% 


87.5±10.4 


(4.26±0.80)xl0" 


-4 


0.187±0.010 


0.170±0.008 


4.31/29 


5x3 


30-60% 


50.3±8.7 


(1.64±0.47)xl0" 


-4 


0.140±0.012 


0.133±0.010 


2.54/29 


6x7 


60-80% 


27.8±2.5 


(1.99±0.43)xl0~ 


-5 


0.0731±0.0002 


0.0719±0.0002 


9.91/29 


5x6 


80-92% 


10.0±1.2 


(1.24±0.30xl0~ 


5 


0.0565±0.0002 


0.0558±0.0001 


8.71/29 


12x12 


STAR 0-5% 


477±13 


(1.48±0.05)xl0" 


'4 


0.140±0.001 


0.132±0.001 


56.6/31 


6x6 


5-10% 


443±15 


(1.08±0.06)xl0" 


-4 


0.127±0.002 


0.122±0.002 


38.0/31 


7x6 


10-20% 


326±17 


(1.02±0.10)xl0~ 


-4 


0.126±0.004 


0.121±0.003 


33.8/31 


6x5 


20-30% 


236±14 


(8.15±1.00)xl0~ 


-5 


0.119±0.004 


0.115±0.004 


30.0/31 


6x6 


30-40% 


169±10 


(7.13±0.09)xl0" 


-4 


0.113±0.004 


0.109±0.004 


25.5/31 


6x6 


40-60% 


109±4 


(4.40±0.22)xl0" 


-5 


0.0961±0.0014 


0.0937±0.0013 


24.5/31 


6x7 


60-80% 


46.0±1.0 


(2.80±0.08)xl0" 


-5 


0.0797±0.0001 


0.0782±0.0000 


23.4/31 


6x7 


PP (nsd) 


4.98±0.15 


(2.87±0.08)xl0- 


-5 


0.0725±0.0000 


0.0711±0.0000 


38.2/28 


20x22 



